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ABSTRACT 
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rn 

^ ! Standard Big Bang Nucleosynthesis at the baryon density determined by the 

microwave anisotropy spectrum predicts an excess of ^Li compared to observa- 
tions by a factor of 4-5. In contrast, BBN predictions for D/H are somewhat 
below (but within 2cr) of the weighted mean of observationally determined val- 
ues from quasar absorption systems. Solutions to the ^Li problem which alter 
the nuclear processes during or subsequent to BBN, often lead to a significant 
increase in the deuterium abundance consistent with the highest values of D/H 
seen in absorption systems. Furthermore, the observed D/H abundances show 
considerable dispersion. Here, we argue that those systems with D/H ~ 4 x 10~^ 
may be more representative of the primordial abundance and as a consequence, 
those systems with lower D/H would necessarily have been subject to local pro- 
cesses of deuterium destruction. This can be accounted for by models of cosmic 
chemical evolution able to destroy in situ Deuterium due to the fragility of this 
isotope. 
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Introduction 



The importanc e of ^Li to cosmology began with the discovery of the Spite plateau 
(ISpite &: Spitelll982l ) . At a time when there was still considerable uncertainty in the baryon 
density, lYang et al.l ( 119841 ) derived an upper limit to the baryon-to-photon ratio, r], effectively 
using the ^Li/D ratio. Remarkably, the ^Li plateau has m aintained a very constant value 



over t he last 30 years and over many obseryatioris, see e.g., (ISpite et al. 



1986 



Hobbs fc Duncan 1987: Rebolo et al. 1988; Hobbs fc Thorburn 



: ISpite &: Spite 



1984 



1991 



Spite &: 



Thorbur: 



3 



Spite 



1994 



Pilachowski et al. 



1993 



1993; 



Molaro et al. 



1995 



Bonifacio fc Molaro 



1997; 



2010 



Rvan et al.lll999uAsplund et al.ll2006l : iBonifacio et al.ll2007l : iHosford et al. 



2009 



2010 



Sbordone et al 



Spite fc Spitell2010l ). While the standard big ban g nucleosynthesis (SBBN) predictions 



for ^Li abundances have also remained relatively stable (IWalker et al 



locco et al 



Cvburt et al.kOOlHCoc et al. 



2002 



2009[ ICoc et Vandoni 



2004 



2010 



mined to unprecedented precision ( iDunkley et al. 



:ICvburt 



2004; 



Coc et al 



1991 



Cuoco et al.ll2004J : ICyburt et al. 



Olive et al. 



2pod; 



2008.: 



20121 ) , the baryon density has been deter- 



2009t iKomatsu et al.l 1201 ll ) from analyses 



of microwave background anisotropies. This has led to a more precise predicti on of the 
^Li abundance that points to a clear discrepancy ( ICyburt et al.l 120081 : lFieldsll201ll ) between 
theory and observation. 



At a baryon-to-photon ratio of ?7 = 6.16 x 10 ^° (iKomatsu et al.ll201ll ). the BBN pre 



diction for ^Li/H is (5.07; 



-0.71\ 



xlO 



-10 



m 



Cvburt et al.l (l2008l ) and 5.24 x 10"^° in lCoc et al 



( 120121 ) with an estimated error bar of 0.5, whi ch is considerably higher than almost all ob- 



servational determinations. The value found in iRvan et al.l (|2000[) was Li/H = (1.23_o'i6) x 



10 ^°. Similarly, the recent analysis of ISbordone et al.l (l2010l ) found ^Li/H = (1.58 ±0.31) x 
10~ ^°. Li observations have a l so be en performed in some globular clusters. For exam- 
ple, IConzalez Hernandez et aP toO^ found (2. 34 ± 0.05) x lO"^" iri NGC 6397, somewhat 
higher than the result for field stars, whereas iMonaco et al.l (120101 ) found a value ^Li/H 
= (1.48 ± 0.41) X 10^^*^ in Omega Centauri that is similar to halo star abundances. 

Resolution of the ^Li problem has involved many different approach es. These range 
from questioning the nuclear reaction rates used in BBN calculations (ICoc et al.l |2004J : 



Angulo et al 



2005 



Cvburt et al 



2004 



Bovd et al.ll2010 ) , or considering additiona l resonance 



19981: 



Pinsonneault et a 



Garcia Perez et al. 



19981 I2OO2I : iRichard. Michaud fc Riched l2005l : iKorn et al. 



reactions ( ICyburt fc Pospelovll2009l : IChakraborty et al.ll201ll : iBroggini et al.ll2012|). The pos- 
sibili t y that depletion plays a role has been discussed at length (IVauclair fc Charbonnel 



20061 : 



20081 ). The temperature scale used in the ^Li abundance det ermination 
has also been considered ( jMelendez fc RamirezI I2004J : iHosford et al.l l2009l I2OIOI ). There is 
also the possibility that the solution of the ^Li problem requires physics beyond the stan- 
dard model. For example, the decay of a massive particle during or after BBN could affect 
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the light element abundances and potentially lower the ^Li abundance (iJedamzikl 12004 



Kawasaki et al 



2005 



Fenget all 12004 : 



Ellis et al. 20051: Jedamzik et 



3 



al. 



2006 



Cvburt et al. 



200f 



Kusakabe et al 



12007 



a 



Cumberbatch et al. 



1200 



Kawasaki et al. 



'ospelovll2009[ICvburt e[a 



20081: 



Pospelov et al. 



20081: 



Jittoh et al 



Jittoh et al. 



20081 



2010l: ljedamzikll2008alJb : 



2010; 



edamzik fc 



Baillv et al. 



20091: 



Pospelov fc Pradlei 



Kawasaki fc Kusakab 



J boiih . 



2009 



Kusakabe et a. 



2010aB : Icvburt et al. 



2008 . 



201 



a 



Another possibility raised recently is that 
of an axion condensate which cools the photon backgro und leading to lower value of y at 



the t ime of BBN relative to that determined by WMAP (JErken et al.ll2012l : iKusakabe et al. 



2OI2I ). More exotic solutioiis involve the possibility of a variatioii in the fundamental con- 
stants JPmitriev et al.ll2004J : ICoc et al.ll2007l : iBerengut et"aDl201o[ ). 



The problem with ^Li should be put into context with respect to the other light ele- 
ments produced in BBN. In particular, there is relatively good agreement between the BBN 
predictions for ^He and D/H and their observational determinations. The helium abundance 
is the most accurately predicted of the prim ordial abundances. A t ?7 = 6.1 6 x 10~^°, the 



helium mass fraction is Yp = 0.2483 ± 0.0002 (ICvburt et all (120081 )): ICoc et all (120121 ) found 



a similar value, Yp = 0.2476 ± 0.0004. On the observational side, the determination of the 



heliu m abundance in extragalactic HII regions is plagued with difficu lties (lOlive fc Ski 



20011 ). Using the Markov Chain-Monte Carlo rn e thods described in lAver et al.l (1201 ll ) and 



Iman 



data compiled in llzotov et al.l (120071 ) . lAver et al.l (120121 ) found a higher value with consider- 
ably larger uncertainties, Yp = 0.2534±0.0083. Given the uncertainty, this value is consistent 
with the BBN prediction. 

The observationally determined deuterium abundance is also in reasona ble agreement 
with it s BBN prediction. There are nine quasar absorption system observations (IBurles fc Tvtler 



1998a 



2003 



b; O'Meara et al. 



O'Meara et al. 



200 



2001 

T 



Pettini fc Bowen 



Pettini et al. 



I2OO8 



2001 



Levshakov et al. 2002; Kirkman et al. 



Fumagalli et al.ll201ll ) with measurable D/H 



and with a weighted mean abundance of3.05±0.22xl0 ^. This shoul d be compared to the 
BBN prediction at the WMAP v a lue o f r] of 2.54 ± 0.17 x 10"^ from ICvburt et all J2008h 
and 2.59 x 10~^ from ICoc et al.l (120121 ) with an estimated error of 0.15. The individual 
measurements of D/H show considerable scatter (a sample variance of 0.62) and it is likely 
that systematic errors dominate the uncertainties. While the agreement is certainly reason- 
able, we do draw attention to the fact that the predicted abundance is somewhat lower than 
the observed mean and this (slight) discrepancy will go in the direction of the possible ^Li 
solution discussed below. 

As might be expected, when one attempts to resolve the ^Li problem by going beyond 
the standard model, there may be consequences for the other light elements. Indeed, the 
destruction of ^Li is often accompanied by the production of deuterium. In the models we 
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will consider below, we assume a relatively late decay of a massive particle. Specifically (to 
be described in more detail in the next section), we consider the decay of a massive gravitino 
(mass 3-5 TeV) with a life-time of order 100-500 s. The decays eject both hadronic and 
electromagnetic energy which breaks up a small fraction of ^He leading to an increase in the 
deuterium abundance and free neutrons whic h destroy the f r eshly pro duced ^Be from BBN . 
The models considered here were developed in lCyburt et al.l (120091 ) . In lCyburt et al.l (120101 ) . 
the particular part of the parameter space most efficient for destroying ^Li was studied and 
we make use of those results here. 

In what follows, we will briefly review the results of post-BBN processing by a late 
decaying massive particle. We will see a tight correlation between the D/H and ^Li/H 
abundances. As a result, a solution to the ^Li problem will invariably lead to an excess of 
D/H. However, deuterium can be easily destroyed in stars. In section 3, we first attempt 
to trace the evolution of these abundances in a model of cosmic chemical evolution based 



on hierarchical clustering (JDaigne et al.ll2004l . |2006| ). As we will see, constraints from the 
cosmic star formation rate will limit the average amount of destruction of both D and Li. 
Furthermore, we will argue that the observed scatter is probably due to the in situ destruction 
of D/H. In section 4, we analyze the deuterium data and in section 5 we discuss new and 
interesting observations of HD/H2 and the possibility that only those absorption systems 
with the highest D/H represent the (post)-primordial value. 



2. Post-BBN processing of the light elements 

While there are many possible models which lead to the post-BBN processing of the 
light elements, we will focus here on a set of supersymmetric models in which there is a 
massive gravitino which deca ys into the lightest supersymrnetric particle (LSP) which is 



also a dark matter candid ate (|Jedamzikll2004l : iKawasaki et al.ll2005l . l2008l : lKohri et al.ll2006 



Cyburt et al.ll2009l . |2010| ). In particular, we consider a constrained version of the minimal 
supersymmetric model which is described by four parameters: a universal gaugino mass, 
mii2'i a universal scalar mass, mo; a universal trilinear term, Aq\ and the ratio of the two 
Higgs expectation values, tan/3. In addition, we must specify the sign of the Higgs mixing 
param eter, /i, which take to be positive, and the gravitino mass (see e.g. (JEUis fc Olive 
2010h l 



As we will see, there is very little dependence on the choice of the specific supersym- 
metric model defined by the parameters r7ii/2,''^0; ^o? and tan/3. However, resulting abun- 
dances of D and Li will be quite sensitive to the the gravitino mass a nd its assumed abun- 



dance. We will focus on a small but varied set of benchmark points (JBattaglia et al.ll2001 
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2004 : Ide Roeck et al.l 120071 ) which are defined by (mi/2,mo, tan/3) = (400,90,10) - point C; 
(300,1615,10) - point E; (460,310,50) - point L; (1840,1400,50) - point M; all with Aq = 0. 
The post-BB N processing of agrayitino decay for each of these benchmark points was studied 
extensively in ICyburt et al.l (120091 . |2010| ) and we make use of those results here. 



In the models of interest, gravitinos decay to lighter supersymmetric states, ultimately 
leading to the lightest neutralino in the final state. For the post-BBN processing of the light 
elements and the destruction of ^Li, decays with the injection of hadrons are most important. 
These occur through gravitino decays to a neutralino + Z, or when kinematically allowed to 
a chargino + W, and gluino + gluon. Three-body decays which include a quark-antiquark 
pair may also be important. Subsequent decays of the gravitino decay products (e.g. Z 
decays) lead to the injection of non-thermal protons and neutrons which may interact with 
the recently produced light elements. 

The dominant primary reaction involving the non-thermal decay products is the photo- 
erosion and spallation of ^He. Since ^He dominates the abundances of the other light elements 
by many orders of magnitude, even a small amount of destruction (which leaves little trace 
on the ^He abundance itself) may induce significant changes in the abundances of the other 
light elements. Generally, the spallation of ^He leads to a net increase in the deuterium 
abundance, as one might expect. The final abundance of ^Li depends sensitively on the 
lifetime of the decaying particle. For relatively long lifetimes (> 10^ s), the ^Li abundance 
is increased. The A=3 products of ^He spallation, tritium and ^He undergo secondary 
interactions ^He(a,7)^Be and T(a,7)^Li, both of which increase the final ^Li abundance. 
For shorter lifet imes, there is a two-step process which leads to the net destruction of ^Li 
(j JedamzikI |2004| ) . First, thermalized and non-thermal neutrons interact with ^Be through 



^Be(n,p)^Li, which is followed by ^Li(p, a)'^He. 



For each of the four benchmark points, ICyburt et al.l (J2010l ) used the observationally 
determined abundances of "^He, D/H and ^Li/H to find the value of the gravitino mass and 
abundance which best fit the data. Indeed, it was found that a significant improvement to 
standard BBN could be achieved with gravitino masses 1713/2 ~ 4 — 5 TeV with abundances 
C3/2 = '"^3/2'^3/2/'"'7 ~ 5 X 10^^^ — 5 X 10~^°. With very little dependence on the specific 
supersymmetric model, the resulting D/H and ^Li/H abundances were approximately 3.2 
xlO~^ and 2.4 xlO^^'^ respectively. The resulting D/H and ^Li abundances for each of the 
four benchmark models is shown in Figure [1] for gravitino masses in the range 2-5 TeV with 
abundances C3/2 = 2.4 x 10"^^ - 2.4 x 10"^ GeV with cuts: ^Li/H < 4 x 10"^° and D/H 
< 10~^. The best fit point in each case, is shown by a (red) star. While there is considerably 
more scatter in benchmark point M, the lower envelope in each of the models is clearly very 
similar. 
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Fig. 1. — The resulting ^Li abundance as a function of the resuhing D/H abundance due to 
the late decays of a massive gravitino. 

While the best fits in each case offer a significant improvement with respect to the 
observed abundance^, from Figure [H one can see that the ^Li abundance can be lowered 
further at the expense of slightly higher D/H. However, the D/H abundance is already 
relatively high at the best fit points with respect to the weighted mean of D/H as determined 
from observations of quasar absorption systems as described in the previous section. We next 
examine the role of cosmic chemical evolution on these abundances and take a critical look 
as to whether the weighted mean is an appropriate measure of the primordial deuterium 



^The x^ based on ^He, D/H, and ^Li/H for standard BBN is 31.7 . In contrast, the y ^ for benchmark 
points C,E,L,M are 5.5, 5.5, 5.4, 7.0 at their respective best fit points (jCvburt et al.ll201Cll ). 
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abundance. 



The effects of cosmic chemical evolution 



3.1. Generalities 



To follow the cosmic chemical evolution of the light element s along with a metallicity 
tracer, we use an analytical model developed first by lDaigne et al.l ( 12004 l2006l ). reproducing 
the cosmic star formation rate in a cosmological context of structure formation. T he model is 
based on the standard Press-Schechter (PS) formalism (jPress fc Schechterlll974l ) to account 
for non linear structures. The rate at which structures accrete mass is determined by a 
Press-Schechter distribution function, fps{M,z). The model tracks baryons 1) within stars 
or their remnants within collapsed structures, 2) in gas within collapsed structures (the 
interstellar medium, ISM), or 3) outside of structures (the intergalactic medium, IGM). The 
model includes mass (baryon) exchange between the IGM and ISM, and between the ISM 
and the stellar component. The age t of the Universe is related to the redshift by 



dt 

dz 



9.78h-' Gyr 



(1) 



1 + Z)x nA + nra(l + z) 



assuming the cosmological parameters of the so-called "concordance model" , with a density of 
matter Q^ = 0.27 and a density of "dark energy" Q\ = 0.73 and taking Hq = 71 km/s/Mpc 
for the Hubble constant {h = 0.71). This allows us to trace all of the quantities we describe 
as a function of redshift. The input stellar data (lifetimes, mass and type of remnant, metal 
yields, and UV flux) are t aken to be dependent on both the mass and the metallicity of 
the star (see iDaigne et al.l (120041 ) for more detail). Once these parameters are specified, we 
can follow many other astrophysical quantities such as the global SFR or the abundances of 
individual elements (Y, D, Li, Fe, O) in the course of the expansion of the universe. This has 
been u sed to investigate sp ecific issues related to early star formation and the reionization 
epoch jRoUinde et al.lbood l 



In this study, we will consider a combination of three distinct modes of star formation: a 
normal mode of Pop II/I, and two additional modes of m assive and intermediate mass (IM) 
Pop III stars (for more detail see IVangioni et all (1201 ll )). Each mode has a specific IMF: 
between 0.1 Mq and 100 Mq for the normal mode of star formation, between 2 Mq and 8 
Mq for the intermediate mass mode and between 36 Mq and 100 Mq for the massive mode 
; the slope of the IMF is taken to be close to the Salpeter value {x = 1.6). 



There are several motivations for including the intermediate mass mode. These stars 



may c orrespond to Pop II I. 2 stars which originate from material polluted by pristine PopIII.l 
stars (JBromm et al.ll2009l ). However, our mass range precludes any type II supernovae (SNe) 
associated with this mode of star formation. In addition, there are theoretical arguments that 
the metal-free IMF predicted from opacity-limited fragmentation theo ry would peak aroun d 
4-10 Mq with steep declines at both larger and smaller masses (JYoshii fc Said Il986l ). 
Primordial CMB regulated-star formation may also lead to the production of a population 



of ea r ly intermediate mass stars at low metallicity ( ISmith et al. 
20ld ISafranek-Shrader et al.lboioh . 



20091 : ISchneider fc Omukai 



There is considerable evidence for an early contribution by IM stars from observations. 
These stars produce very little in the way of heavy elements (oxygen and above), but produce 
significant amounts of carbon and/or nitrogen and above all helium. Ev idence exists that 
the number of carbon-enhanced stars increases at low iron abundances ( JRossi et al.lll999[ ) 
necessitating a Po p III source of carbon, possibly in the asymptotic giant bra nch (AGB) 
phase of IM stars (JFujimoto et al.l l2000l : lAoki et al.ll2002l : iLucatello et al.l 120051 ) indicating 
possibly an IMF peaked at 4 - 10 M© (JAbia et al.l 1200 ll ). In addition, the presence of s- 
process eleme nts, particularly Pb at very low met allicity also points to AGB enrichment 
very early on (jAoki et al.l l2001t ISivarani et al.l |2004J ) . 



Isotopic studies of Mg also show a need for an early generation o f IM stars. While c ore 
collapse supernova e produce almo s t excl usively ^^Mg, observations of lYong et al.l ( l2003al Jbl): 



Alibes et al.l ( 120011 ): iFenner et al.l (l2003l ) show enhancements (relative to predictions based 
on standard chemical evolution models) in both ^^'^^Mg. Finally, IM stars may in part be 
responsible for the s omewhat high ^He abundances seen in low metallicity dwarf galaxies 
(IVangioni et all 120111 ). Our motivation here for considering an early population of IM stars 
is clear: they offer the pos sibility for destroy ing significant amounts of Deuterium without 
producing heavy elements (JFields et alll200ll ). 



We fit the SFR history of Pop II/I stars to the data compiled in iHopkins &: Beacom 



Bouwens et al. 



( 2006 ) (fro m z = to 5), and to the recent measurements at high redshift by 

(2010) and lGonzalez et al.l ( 120101 ). This is modeled by the expression suggested by lSpringel and Hernquist 

(l2003l ) : 

a — + oexp(a [z — Zm)) 

where i> and Zm are the amplitude (astration rate) and the redshift of maximum SFR re- 
spectively; and b and b — a are related to the slopes of the curve at low and high redshifts 
respectively. The normal mode is fit using: z/n/i = 0.3 Mq yr~^ Mpc~^, -Zmii/i = 2.6, 



flii/i = 1-9 and 6n/i = l-l- The SFR of this mode peaks at z 
strong constraints on the Pop II/I SFR. 



3. These observations place 
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Fig. 2. — heit panel: Cosmic star formation rate (SFR) as a function of redshift. The data 
(solid black points and solid error bars) are taken from Hopkins fc Beacom ( 20061). Dashe d 



black error bars correspond to results by iBouwens et al 



( l2010l ) and lConzalez et aP toid ). 



The blue solid line represents the normal SFR mode with a Salpeter IMF and a mass range: 
0.1 < M/Mq <100. The dotted red line represents the massive Pop III stellar mode with 
a mass range: 36 < M/Mq < 100. The dashed black curve represents the intermediate 
mass SFR mode, 2< M/Mq < 8. Right panel: the optical depth deduced from WMAP 
observations is presented as a fu nction of the redshift . The red range corresponds to the 
observed results from WMAP7 ( jKomatsu et al.l ( 1201 ll )). The red line corresponds to the 
model including all SFR modes. 



As noted above, in addition to the normal mode, we add two modes for Pop III stars. 
The IM (massive) SFR parameters are: z/ma = 0.4(0.7) Mq yr~i Mpc~^ Zmiiia = 2.8(10.), 
fliiia = 1.9(4) and fema = 1.1(3.3). Parameters are chosen to maximize the potential for D/H 
astration (without the production of heavy elements) while remaining conco rdant with the 
observ ed S FR at high re d shift. For a detailed description of the model see iRoUinde et al. 
( 120091 ) and iDaigne et al.l (120061 ). A recent study show s that the first stars have feedback 
limited masses of around 40 M© (JHosokawa et al.ll201ll ). appropriate for core-collapse super 



novae. 



In Figure HI we show the adopted SFR for each of the three modes considered. As one 
can see, the normal mode and the intermediate mass SFR are constrained to fit the obser- 
vations plotted in the figure. As the data extend only np to z ^ 8, there is little constraint 
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for z > 8. Indeed, data at these redshifts are highly uncertain due to unknown systematics 
invol ving, among other effects, dust correction and adopted rest-frame UV luminosity func- 
tion (JLabbe et al.ll2010f ). In this context, the massive mode which can dominate at 2 > 10 is 
constrained to fit the optical depth deduced from WMAP observations (shown in the right 
panel of Figure [2]). The global SFR used in this study is the sum of these three modes. 



3.2. The cosmic evolution of D and Li 

In the cosmological context described above, it is possible to track the evolution of both 
light elements, D and ^Li. Following the analysis of Section 2 corresponding to the post-BBN 
processing of the light elements, we consider the resulting D and ^Li abundances due to the 
late decays of a massive gravitino (see Figure 1). As noted earlier, a general consequence 
of lowering ^Li is higher D/H. So to fit the Spite plateau we choose three representative 
values for the ^Li abundance corresponding to the values quoted in the introduction: ^Li/H 
= 1.23 X 10~^°, 1.58 X 10"^°, and 2.34 x 10"^°. These abundances correspond to D/H = 
4.4 X 10-^ 3.9 X 10~^ and 3.3 x 10"^ respectively. 







BBN value 



Mean value, Sbordone et al (ZOT?!)'' 



.J.,::... 



-2.5 

[Fe/H] 



Fig. 3. — D/H abundance as a function of redshift (left) and A(Li) = log{Li/H) + 12 as a 
function of [Fe/H] (right ). Deuterium data co me from observations described in the text. 
Lithium data come from (jSbordone et al.ll2010l ). The red solid (black short dashed and large 
dashed) lines correspond to initial post-BBN Li/H, D/H values: 1.23 x 10"^°, 4.4 x 10"^, 
(1.58 X 10-1° , 3.9 X 10-^ and 2.34 x 10"^°, 3.3 x 10"^) respectively. 
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In Figure [3] we show the cosmic evolution of hthium as a function of metalhcity and the 
corresponding D/H evolution as a function of redshift in the three selected cases. At the 
metallicities of interest, there is very little evolution in the ^Li abundance. The high redshift 
Deuterium data come fror n observatiori s discu ssed in the next section. The local Deuterium 
data (z = 0) i s take n from lLinsky et al.l ( 120061 ) in the Local Bubble and the higher value from 



Savage et al.l ( 120071 ) in the warm neutral medium of the Galactic halo. D destruction begins 
at redshift 3, corresponding to the peak of the cosmic SFR (see figure 2, left). The global 
evolution of D/H fits the observations reasonably well, given the observational uncertainties 
(see next Section for a discussion). 

It is important to note the high dispersion of the observed values at 2; ~ 3 (correspond- 
ing to t = 3—4 Gyr) in the range 2.4 — 4. x 10^^. This could be a consequence of different 
star formation histories in the galaxies associated with the DLAs. Note that the individ- 
ual chemical evolution histories are not probed by our model of cosmic chemical evolution 
which only yields the average abundances of the elements in structures and in the IGM. 
However, D is a very fragile isotope (it is destroyed in stars at T = 10^ K) and its destruc- 
tio n rate is highly dependent o n the ratio a = gas mass/total mass of the galaxy as shown 



m 



Vangioni fc Audouzd (119881 ) (see their figure 3); a low ratio corresponding to a low D 



abundance. As a consequence, D/H destruction factors can range between 2 to 10. Indeed, 
Fields et all (I2OOII ) showed that D can be efficiently destroyed without overproducing metals 
in a standard chemical evolutionary model with an early population of intermediate mass 
stars (see their figure 1). It is reasonable to think that these pristine structures are in dif- 
ferent evolution stages and consequently have different gas masses and processing histories . 
Finally, as D can only be destroyed in the course of chemical evolution ( lEpstein et al.lll976l ). 
it is reasonable to relate the highest D/H observation to the primordial post-BBN abundance. 
In this context, we are able to reconcile the BBN D abundance and high D/H observational 
values at high redshift together with theoretical and observational Li abundances. 



4. Analysis of the D/H observations 



Most of the measurements available in the 



(I2OO8I ) and more recently in iFumagalli et al.l (120111). It was concluded, however, that the 



iterat ure have been gathered by lPettini et al. 



measurements by iLevshakov et al.l (120021 ) and ICrighton et al.l (120041 ) are not relevant because 
not all of the D I components in those observations are resolved. This effect would imply 
systematically low D/H if one has included velocity components in H I that are not dete cted 
in D I because the D I is hidden by other H I absorption features ( JKirkman et al.ll2003l ). 
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Levshakov et al.l (120021 ) convincingly showed that D I is detected at z = 3.03 toward 
Q 0347—3819 in the transitions Ly-8, 10 and 12 (see their Fig. 12). They fit the D I absorp- 
tion with a profile that is consistent with the N I and O I profiles. Although the latter is 
saturated, the former is well defined in N lAA953.4,964.0 (see their Fig. 10). The high D l/H I 
value derives from the presence of a narrow component with Doppler parameter as small as 
~3 km s~^. The presence of this narrow component is well demonstrated however by the de- 
tection of H2 in a single component. Thirty one H2 transitions are used to derive the Doppler 
parameter from the curve of growth (see their Fig. 3). In addition the H2 excitation diagram 
indicates a temperature of ~825 K which is consistent, within errors, with the low value 
of the Doppler parameter. The temperature found here is also consistent w ith the findings 



that the molecular fraction is small in most DLAs because the gas is warm (IPetitiean et al. 



20061) . Although errors may have been underestimated (see Appendix of iKirkman et al. 
( 120031 ). there is no objective reason to reject this me asurement, D/H =( 3.75x0. 25) xlO~^ 
for log A^(H l) = 20.56±0.05 and Zn/H = 0.98±0.09 dhedoux et al] l2003h. Thus, we have 



included this value in our study. Note that iD'Odorico et al.l (I2OOII ) report a lower value in 
this system: D/H = 2.24±0.67xl0~^, but their b values are certainly too large and their 
value should be considered as a lower limit. 



Crighton et al.l (120041 ) reported D/H = 1.6+!]:f xlO^^ towards PKS 1937-1009 in a Ly- 
man Limit System (LLS) with log iV(H l) = 18.25±0.02 and for a metalhcity of [Si/H] = -2.0± 
0.5 dX z = 3.256. In that case only Lyman-a and Lyman-/? lines are used and i ndeed, a model 



with l ower H I and higher D I column densities can be accommodated (see lO'Meara et al. 



( 120061 )). As a result, we have not included this reported measurement in our analysis. 



In addition, i t is su rprising that iPettini et al.l ( l2008l ) do not consider the measurement by 
Levshakov et al.l (120021 ) when they retain the measurement at 2; = 2.06 towards Q 2206—199 
which is obtained from a low SNR intermediate re s olutio n {R ~ 15000) STIS-HST spectrum 
(see discussion in the Appendix of iKirkman et al.l (120031 )). Furthermore, the h value and the 
decomposition of the profile are uncertain in t his observation and t he D/H ratio has been 
claimed to be possibly as high as 2.9x10"^ (see (IKirkman et al.ll2003l )). While we retain this 
object in our average, we comment as well on the effect of excluding it. 



Very recently, iFumagalli et al.l (l201l[ ) measure log (D/H) = — 4.69±0.13 in a LLS 
(log N{Y{ l) = 17.95) aX z = 3.522 towards SDSS J113418.96+574204.6. They argue that 
the gas is pristine with metalhcity <10~^'^ solar. However, at such a low neutral hydrogen 
column density, the gas is most probably partially ionized and a correction for ionization 
effects has to be applied. To reach their conclusion, they assume a value for the ionization 
parameter (log U > —3) and argue that a gas density larger than njj > 10"^ cm~^ is unusual 
in LLSs. We believe that this assumption is unjustified and a density of njj ~ 10^^ cm"'^ 
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is perfectly acceptable (see e.g. iPetitjean et al.l (Il992[ )). Therefore a robust upper limit on 
the metallicity is not less than 10~^ solar. This abundance is typical of w hat is found for 



Darnped Lyman a Syste ms (DLAs) with the smallest observed metallicities (jPenprase et al. 
20 111 : ICooke et al.l |20 111 ). Therefore, the gas may not be pristine and the low value of the 



D/H ratio could be a consequence of deuterium destruction by star formation activity. 



Consequently, we hav e included all deuterium measurements from lPettini et al.l (120081) 



plus t he measurements by iLevshakov et al.l (120021 ) and the recent report in iFumagalli et al. 



(120111 ). The weighted mean of these nine measurements is D/H = (3.05 ± 0.22) x 10~^, 
where a scale factor of S* = 2.0 (= a/xVS) is applied to the the error in the mean to 
partially account for the large dispersion. The sample variance of would imply a much larger 
uncertainty of 0.62 xlO~^. We n ote that this va l ue is in modest disaQreement with the SBBN 



value of 2.54 ± 0.17 x 10^^ from jCYburt et al.l J2OO8I) o r 2.59 ± 0.15 x 10"^ from ICoc et al 



(120121). Excluding the data from ILevshakov et al.l (|2002|), we find D/H = (2.80 ±0.20) x 10"^ 
{S = 1.64) and a sample va r iance of (0.52 xlO~^) for eight objects. Neglecting instead the 
data from IPettini fc BowenI J200lh . we find D/H = (3.11 ± 0.21) x 10"^ {S = 1.83) and a 
sample variance of (0.55 xlO~^), also for eight objects. 

It has already been noted by several authors that the dispersion in the reported D/H- 
values is much larger than what is expected from individual errors. This could be a conse- 
quence of the errors having been underestimated. The presu mption of a uii i que v alue for 
D/H, however, is not supported by the observations (see also llvanchik et al.l (120 lOh ). Fur- 
thermore, if the post-BBN value for D/H were at or near the upper end of the existing data 
as in Figure m the dispersion seen in the data could be explained by the in situ destruction of 
D/H through chemical evolution. If this is indeed a local effect, the extent of the dispersion 
is not surprising. 



5. Discussion 



As mentioned in the Introduction, there are many potential solutions to the SBBN 
^Li problem. It is fair to say that all of them require some additional input beyond what 
we normal term as standa rd. This may in volve additional turbulent features in diffusive 



models of stellar evolution (iKorn et al.ll2006f) or new resonant reactions (ICyburt fc Pospelov 



2OO9I : IChakraborty et al.ll201ll : iBroggini et al.ll2012l ). However, several of the possible solu- 



tion come with the pr ice of an increased D/H abundan ce: photon cooling (lErken et al 



Kusakabe et al 



2005 



20121) . variable fundamental constants (iDmitriev et al. 



Berengut et al.ll2010l) and particle decays during or after BBN (iJedamzik 



Feng et al.ll2004l : lEUis et al.ll2005l : Ijedamzik et al.ll2006l : ICvburt et al.ll2006 ; 



2004 



Coc et al 



2004 



2012 



2007 



Kawasaki et al 



Kusakabe et al. 
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2007: Cumberbatch et al.ll2007: Kawasaki et al 



Jedamzik fc Pospelovlbood: ICvburt et alJboOQi 



2008: Pospelov et al 



Kusakabe et al.ll2008 



2008: Jittoh et al. 2008 



2010 



Jedamzik 



2008a 



Bail 



V et alJbood : IPospelov fc PradleJboiOaE I: ICvburt et allboiot Ijittoh et al.lboiol : 



Kawasaki fc Kusakab 



20 111 ) all tend to increase D/H while yielding a ^Li abundance which matches t he low metal- 



l icity Pop II abundance determinations. Here, we have used the results of ICyburt et al. 



( 120101 ) for the specific correlations between the post-BBN D/H and ^Li/H abundances as 
displayed in Figure [H 

Since the D/H abundance is expected to decrease monotonically over time, a high value 
of D/H from BBN is not necessarily problematic. First of all, we emphasize that SBBN 
predicts a value of D/H which is about 2cr below the weighted mean of observational deter- 
minations. Second, while the post-BBN abundance of D/H may exceed the observational 
value, we expect some destruction of deuterium due to chemical evolution. Indeed in a model 
based on hierarchical clustering, we have seen that there is some modest astration of D/H 
while leaving ^Li virtually unperturbed as seen in Figure [31 However, because these models 
yield average abundances, they can not account for the observed dispersion seen in the data. 
As noted above, this may be due to either under-estimated errors or the in situ destruction 
of D/H. 

Before concluding, we note that the HD/H2 ratio is an interesting alternative for D/H 
investigations in DLAs. In the diffuse ISM, with physical conditions very similar to those 
of DLAs, the formation of HD occurs via the reaction: H2 + D~*" — )■ HD + H"*" while its 
destruction is due to photodissociation. Because of the low abundance of deuterium, the 
transition between atomic deuterium and HD takes place deeper in a cloud than the transition 
between atomic and molecular hydrogen. Moreover, while the H2 column density is usually 
very well constrained by the presence of numerous transitions, only a few transitions are seen 
for HD, making these measur ements quite uncertain . Therefore, although the ratio should 
be interpreted with care (e.g. iTumlinson et al.l (120101)) , iV(H D)/2A^(H2) in the cloud should 
give a lower limit on the D/H ratio (see lLacour. et al.l ( 120051 )). 



Molecular hydrogen is found in less than 10% of DLAs ( Ledoux et al.ll2003l : lNoterdaeme et al. 



2008al ) and to date, only five detections of HD at high redshift have been reported. Interest- 
ingly, some high D/H values are observed despite significant metal enrichment: A^(HD)/2A^(H2) 
= 1.5xl0-^ 3.6xl0-^ 7.9xl0-^ 1.6x10-^ and 0.95 xlO"'^ at z^^^ = 2.42, 2 .33, 2.10, 2.63, 



and 2.69 towards, J 1439+1117 (JNoterdaeme et a 



2OIOI). J 2123-0500, FJ 812+32 (iTumhnson et al. 



j2008aD, Q 1232+082 (llvanchik et al. 

2OI0I) and SDSSJ 123714.60+064759.5 
( iNoterdaeme et al.ll2008bl ) for metallicities relative to solar [X/H] = +0.16, —1.43, +0.36, 
—0.48 and +0.34 respectively. 

It is apparent that these measurements, although uncertain, are all consistent with the 
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cosmological D/H ratio and, m ore importantly, a factor of ten above what is measured in 



our Galaxy fiLacour. et al.ll2005l ). This finding is somewhat puzzling, in light of the diversity 
in the other properties of these systems and in particular their high metallicity. Infall of 
primo rdial gas can help m a intain such high D/H ratios in the course of galactic evolution 



(e.g., iProdanovic fc Fields! ( l2008l )). However, such process may have difficulties alone to 
explain the observations. Consequently, these HD/H2 observations seem to indicate that 
D/H values at hight redshift are somewhat high; indicating again that a high primordial 
deuterium value is privileged. 

In conclusion, we have argued that although many models which are capable of reducing 
the high BBN abundance of ^Li to the level seen in Pop II halo stars come with a price of an 
increase in the D/H abundance, this may well be a feature and not a bug. The weighed mean 
of the D/H abundance data (modestly) exceeds the SBBN value, and carries considerable 
dispersion. If the upper envelope of this data reflects the post-BBN abundance of D/H, 
then these models which attempt to resolve the ^Li problem are preferred. The dispersion 
seen (if real) must then be explained by the astration of D/H in local processes. We note in 
closing that non-standard model explanations of the primordial ^Li abundance generically 
predict ^Li in excess of the standard model prediction as well as an increase in the deuterium 
abundance. Confirmation of ^Li in extremely metal-poor halo stars would provide significant 
support for models in which nuclear processes during or subsequent to BBN resolve the ^Li 
problem. 

We would like to thank Feng Luo for help with the data leading to Figure 1. The work 
of KAO was supported in part by DOE grant DE-FG02-94ER-40823 at the University of 
Minnesota. This work was also supported by the PICS CNRS/USA and the french ANR 
VACOUL. 
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